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Spectrum and wave function of gapless edge modes are derived analytically for a tight-binding 
model of topological insulators on square lattice. Particular attention is paid to dependence on 
edge geometries such as the straight (1,0) and zigzag (1,1) edges in the thermodynamic limit. 
The key technique is to identify operators that combine to annihilate the edge state in the 
effective one-dimensional (ID) model with momentum along the edge. In the (1,0) edge, the 
edge mode is present either around the center of ID Brillouin zone or its boundary, depending 
on location of the bulk excitation gap. In the (1,1) edge, the edge mode is always present both 
at the center and near the boundary. Depending on system parameters, however, the mode is 
absent in the middle of the Brillouin zone. In this case the binding energy of the edge mode 
near the boundary is extremely small; about 10~ 3 of the overall energy scale. Origin of this 
minute energy scale is discussed. 

KEYWORDS: zigzag lattice, effective Hamiltonian, flat edge spectrum, reentrant edge mode 
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1. Introduction 

Recently, a new class of topological insulator (TI) , also 
referred to as quantum spin Hall (QSH) insulator, has 
been attracting much interest in both theory 1-7 ) and 
experiments. 8 ~ 10 ) Different from conventional insulators, 
QSH insulators have topologically protected helical edge 
states with the spectrum lying in the bulk insulating 
gap. 1,12,13 ' Two-dimensinal (2D) topological insulator 
has been realized in HgTc/CdTe quantum wells 8 ) fol- 
lowing theoretical suggestion of Bernevig, Hughes and 
Zhang 6 ) who proposed a useful model, hereafter referred 
to as the BHZ Model. Helical edge state has been further 
studied using the resultant continuum model. 14,15 ' The 
BHZ model is not realistic away from the zone center, 
since it is based on the k ■ p perturbation theory and 
envelope-function approximation. However, a rcgulariza- 
tion of the model using the tight-binding scheme has the 
simple structure for the whole Brillouin zone (BZ), and 
is suitable for studying the general property of QSH sys- 
tems. Note that characterization of the topological prop- 
erty of the system requires information of wave functions 
over the whole BZ. 

In this paper, we derive the spectrum of edge modes 
analytically taking the tight-binding version of the BHZ 
model for the square lattice. In contrast to previous study 
of edge modes, 11 ' 14 ' we take not only the straight (1,0) 
edge but also the zigzag (1,1) edge, and make detailed 
comparison of respective edge modes. This comparison 
is partially inspired by the remarkable difference be- 
tween the zigzag and armchair edges in graphene. 16, 17 ) 
In contrast to the latter, however, the edges in the 
present model should not be taken as representing the 
HgTe/CdTe quantum wells since the short-distance be- 
havior of the model is not realistic. Nevertheless, the 
intuition gained by the exact solution of the simplified 
model should provide useful information for understand- 
ing more complicated systems. 



In a separate paper, 18 ) we have already derived the 
spectrum of both straight and zigzag edges relying par- 
tially on numerical method. In particular, we have found 
for the zigzag edge a reentrant mode with tiny binding 
energy Since the relevant energy is so small as compared 
with the overall energy, it is desirable to characterize its 
nature analytically. Then the origin of tiny binding en- 
ergy should be clarified. In the present paper, we de- 
rive fully analytic expressions of not only the spectrum 
for both edges, but also momentum regions allowed for 
the modes. Namely we derive critical momentum where 
the edge mode merges with bulk excitations. Our ana- 
lytical method is a systematic generalization of previous 
ones 11, 19 ) so that the spectrum can be obtained for gen- 
eral direction of the edge. 

This paper is organized as follows: In §2, we review the 
lattice version of the BHZ model with nearest-neighbor 
transfer, paying attention to its symmetry property in 
the 2D BZ. Sections 3 and 4 arc devoted to analytic 
derivation of edge modes for straight and zigzag edges, 
respectively. We use a systematic method to derive the 
spectrum in terms of the "annihilator" . Existent regions 
of edge modes are derived for both straight and zigzag 
edges. In the zigzag case, we find a novel reentrant edge 
mode with a tiny binding energy below the bulk spec- 
trum. Finally, we summarize the results and discuss their 
implication in Sec. 5. 

2. Model with particle-hole symmetry 

We consider the BHZ model given by the following 
4x4 matrix: 



H(k) 



(1) 



h(k) 
h*(-k) 

where k = {k x , k y ) is a 2D crystal momentum, measured 
from T-point. The lower-right block h*(~k) is a 2 x 2 
matrix, and is deduced from the upper-left block h(k) by 
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time reversal transformation. h(k) is parametrized as, 

dz d x zdy 



h{k) = d(k) ■ er 



where di(k) are given by 

d x (k) = Asmk Xl d y (k) 
d z {k) = A - 25(2 - cosfc 2 



A sin k y 
- cos k y ), 



(2) 

(3) 
(4) 



with the lattice constant a set to unity. ^ We only con- 
sider the case A > 0, B > and A > in this paper. 
The same signs of B and A are necessary for topological 
insulator. The solution for A < is trivially obtained 
from that for H(—k). Each row and column of eq. (2) 
represent spin-orbital states associated with the s-type 
Tq and the p-type Tg orbitals of the 3D band structure 
of HgTe and CdTe. Other parameters which appear in 
Ref.6, i.e., C and D have been set to zero. As a result, 
the spectrum has the particle-hole symmetry that sim- 
plifies the analysis. 

The 2 x 2-matrix Hamiltonian h(k) is equivalent to 
the following tight-binding Hamiltonian on the square 
lattice: 



i, J 



+ ( c\ jt x ci+x, j + c\ ji v c ItJ+ i 



h.c. 



(5) 



where we have introduced for each site (I, J) the two 
component field: 



with orbitals s,p. The tight-binding parameters are given 

by 



£ = (A-4B)a z , 
A 



t, 



1 1, 



A 



Bg 7 _ 



Ba 7 . 



(7) 
(S) 

(9) 



each of which is a 2 x 2 matrix. The down spin part h± 
corresponding to the lower-right block of eq.(l) is ob- 
tained by replacing ti with its complex conjugation t*. 

We consider the BHZ model over the whole BZ of the 
square lattice. Bulk energy E^(k) is written as 

sin 2 k„) 



E b (k) = ± {A 2 (sin 2 k x 

+ [A — 4B + 2B (cos k x 



cos fc a )] 2 | 



1/2 



(10) 



which is symmetric with respect to positive (conduction) 
and negative (valence) energy bands, being the signature 
of the particle-hole symmetry. As one varies mass param- 
eter A, there appear four gap closing points, namely at 
r = (0, 0), X = (tt, 0), X' = (0, tt) and M = (tt, ti). Note 
that these points are invariant against time-reversal op- 
eration. Gap closing at T occurs when A = 0, whereas 
the gap closing at X and X' occurs simultaneously when 
A = 4B, at M when A = 8B. 
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Fig. 1. Straight edge lattice ribbon with two boundaries in (1,0) 
direction. 



3. Straight edge 

3.1 Effective one-dimensional form 

Let us first consider the geometry (see Fig. 1), in which 
electrons are confined to N r rows in a strip between y = 1 
and y = N r , i.e., the two edges are along x-axis. Trans- 
lational invariancc along cc-axis allows for constructing a 
ID Bloch state with a crystal momentum k = k x : 

\k,J) = -±=Ye* kl \I,J), (11) 



'N, 



i 



where /, J represents lattice site, N c is the number of 



(6) sites along the x-axis, and \I, J) 



-i, J 



|0). In order to 



introduce the edges, it is convenient to rewrite eq. (5) 
in form of a hopping Hamiltonian between neighboring 
rows. In terms of the two-component creation and anni- 
hilation operators Cj(fc) and cj(k) associated with Bloch 
state (11), one can rewrite eq. (5) as 



h\ = ^ft-oi(fc), 

fc 

h 01 (k) = Y^c\{k)S{k)cj{k) 
j 

+ Y1 [cj(k)i y c J+1 (k) + h.c. 



(12) 



where £ (k) is given by 

£ (fc) = Asmka x 



(Ab + 2cosfc) o~ Z: 



(13) 



(14) 



Here and in the following, we use the notation Ab = 
A — AB, and take the energy unit so that B = 1. The 
corresponding Schrodingcr equation is given by 

£ (k)*j + ty^j-! + *J * J+1 = £ t (fc)*j (15) 

where Vfj is the two-component amplitude with row in- 
dex J. The straight edges along the J = 1 row and 
J = N r row can be implemented by open boundary con- 
dition * = 3>jv r+1 = 0. 

It has been shown 14 ) in the continuum approximation 
of eq.(15) that there appear modes, which are localized 
on both edges of the system, with the particle-hole sym- 
metric spectrum and a minimum gap at k = 0. The origin 
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of the gap is the overlap of wave functions on different 
edges. As N r tends to infinity, the gap disappears since 
the overlap becomes zero. In this limit one can identify 
edge modes localized on either edge. In the thermody- 
namic limit, the spectrum E(k) of an edge mode near 
J = 1 becomes an odd function of k, and the other edge 
mode has the spectrum —E(k). The continuum approx- 
imation of Rcf.14 has a limited validity, but emergence 
of the zero mode at k = in the thermodynamic limit 
N r — > oo is guaranteed by the time-reversal and particle- 
hole symmetries. 11,19 ' In the following, we deal with the 
thermodynamic limit. 

3.2 Separation into Hermitian and annihilating parts 

Our strategy to obtain the edge spectrum is best illus- 
trated in the straight edge. Although the spectrum in this 
case has already been obtained in the literature, 11 ' 19 ' we 
here present our way of derivation that can be extended 
to the zigzag edge. We try an edge state solution with 
property: = A*,/ = A J+1 * with |A| < l, 19,20 ' and 

derive the vector Then eq.(15) can be written in the 
following form: 



£(k) + At], + }_i y 



y = Poi(\,k)y = E^(k)y, (16) 



where Poi(A, k) can be rearranged as 
P i(A, k) = 7-er, 

with components 

7 X = A sin k 
A 



Jv=i^ (A -A- 1 ) 



7, = (A + A : ) + A B + 2 cos k 



(17) 

(18) 
(19) 
(20) 

For a general complex vector 7, the eigenvalues of P01 
are also complex. However, we obtain real energy E^{k) 
in eq.(16) if one of the following conditions is met: 

(a) all components of 7 are real (including zero) ; 

(b) the nonzero complex components combine to give 
zero when acting on the edge state. Such combination of 
operators is referred to as annihilator. 

The condition (a) is not relevant here, because we would 
then have two eigenvalues ±E-^(k) and corresponding two 
eigenf unctions for a given k. Actually we should have 
only one edge mode for h^(k). Hence we have to accept 
the condition (b), and separate P01 into a Hermitian part 
that gives the spectrum by diagonalization, and the rest 
that makes up the annihilator. 

From eq.(19) and (20), coefficients 7 y and j z are both 
real only if |A| = 1. On the other hand, due to time- 
reversal symmetry and particle- hole symmetry of the sys- 
tem, the eigenvalue E^(k) is an odd function of k in the 
thermodynamic limit. 11,14 ' Therefore, a y and o z must 
both belong to the annihilator, and the only component 
to be diagonalized is a x . The annihilator corresponds to 
either a y + ia z or a y — ia z , depending on the eigenvalue 
±1 of a x . 



Accordingly, we decompose Poi(A, k) as 

Poi(X,k)=H 01 + F 01 , (21) 
Hqi = Asmkcr x , (22) 



FqI = -JyCTy + ~f z O z 



(23) 



where Fqi should form the annihilator. Namely, we im- 
pose the relation 7^ = si^ y , according to the eigenvalue 
s = ±1 of a x . The relation is equivalent to 

(l + ^)A + (l-^)A- 1 +m fc =0, (24) 

where we have introduced the notation: 

m fc = A B +2cosfc. (25) 
By solving eq.(24) with A ^ 2, we obtain 



A±(s) 

with the relation 



-mfc ± yj rrij, + A 2 
2 + sA 



A ± (-,s) = 1/A T ( S ). 
Provided the eigenvalue equation 

i/oi* = E f (k)ty. 



(26) 



(27) 



(28) 



is satisfied, then equality FoiVP = follows with eq.(24). 
In this way the edge spectrum is simply derived as 



E^(k) = s A sink. 
with the eigenstate written as 

T 

.s, 



(29) 



(30) 



Note that only one of s = ±1 is relevant, as discussed 
below. 

The boundary condition ^7=0 = requires the edge 
state to have the form 



(31) 



apart from the normalization factor. Since the wave func- 
tion should decay as J increases, the edge state can only 
be realized with |A±(s)| < 1. Because of the relation 
eq.(27), there is at most one s for given k that describes 
the edge mode with the spectrum eq.(29). Hence there is 
only one edge state per spin and per momentum. 

In a similar manner, we obtain another edge mode for 
ft.4,. The corresponding energy E^(k) is given by 



E^k) - E t (-k) = -E t (k), 



(32) 



where the first equality corresponds to the time-reversal 
symmetry. The particle-hole symmetry connecting the 
rightmost and leftmost sides with the same k involves 
different spins. Note that the Kramers pair has the same 
A±(s) for given k and — k. In certain range of k 7 however, 
the edge modes do not exist. This problem is studied in 
detail in the next section. 

3.3 Allowed momentum range for edge modes 

A pair of gapless edge modes per edge are always 
present in the topological insulator phase (TI) that ap- 
pears for < A < 8, i.e., |A B | < 4 in the BHZ model. 
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Fig. 2. ID energy bands in the straight edge: TI-1 (top) with 
A = 1.2 and TI-2 (bottom) with A = 4.5. In both cases, we 
set A = 1. The edge modes have the spectrum ±7lsinfc in the 
existent range of k. 



However, the number of zero points of helical edge states 
should be an odd multiple of two, including the degen- 
eracy, in the ID BZ. 12 ' If the edge modes with the spec- 
trum E = ± A sink were present for the whole BZ, the 
zero points amount to four (an even multiple of two) 
which violates the topological stability. 21 ) Hence the edge 
modes must merge into bulk excitations at finite K m , and 
degenerate zero points occur either at fc = or k = tt, 
but not at both. Note that fc = and k = tt are two time 
reversal invariant momenta in ID BZ. 

With As — 0, the energy gap closes at X points (w, 0) 
and (0, tt) in the 2D BZ, as seen from eq.(10). Let us 
classify the case of —4 < Ab < as TI-1, and the case of 
< Ab < 4 as TI-2. Figure 2 shows the spectrum of edge 
modes in the TI-1 (upper panel) and TI-2 (lower panel). 
In the TI-1, the edge modes intersect at fc = 0, whereas in 
TI-2 they meet at fc = tt. Also shown is approximate bulk 
spectrum in the system that is derived from E^k, k y ) by 
fixing k y to 2Trn/N r with n = 1, . . . , N r . Each curve for 
the bulk spectrum corresponds to integer n. 

We now derive the pair ±fc m of momentum where the 
edge mode merges with bulk excitations. In the follow- 
ing, we always assume k > for simplicity. The merging 
occurs when the larger of |A±(s)| becomes unity. Note 
that such A±(s) is real, since otherwise |A±(s)| is inde- 
pendent of rrifc, and hence of k. For complex A±(s) we 
obtain from eq.(26), 

|A ± (.s)| 2 = §-^, (33) 

which is less than unity only with s = 1. In this section 
we deal with the case < A < 2 that allows complex 



A±(s). Then, edge modes are assured to be present for 
such k with complex A±(l), which is simply written as 
A± hereafter. Edge modes in the case of A > 2 will be 
discussed in §5. 

TI-1 regime 

In this case gapless points are present at k = 0. Let us 
assume mu = Ab + 2 cos k < when merging occurs at 
k = k m i, Then we obtain |A + | > |A_| for real A±, and 
the condition for merging is reduced to 



A+ = 



-rrtk + \J ml + A 2 



2 + A 

from eq.(26). Then we obtain 

cosfc m i = 1 - A/2, 



= 1. 



(34) 



(35) 



which justifies the assumption mk < 0. One can check 
that there is no solution if we assume m k > 0. We note 
that eq.(35) has already been obtained by Konig et al. 11 ' 
The condition for merging is also to have the same 
energy as the lowest bulk excitation for given k. The 
minimum of Et(k, k y ) may occur either at k y = or k y = 
iarccos [2rrik/(A 2 — 4)] depending on the value of fc. If 
the threshold of bulk excitations with fc = k m i occurs at 
k y = 0, merging momentum k m \ is simply obtained from 
eq.(10) as 



A B + 2 + 2cosfc = m k + 2 = 0, 



(36) 



which is consistent with eq.(35). Namely, the condition 
A + = 1 is equivalent to having the same energy for edge 
mode and for the minimum of bulk excitations. Further- 
more, it is easily seen that the group velocity A cos k m \ at 
merging point is common to both edge mode and the low- 
est bulk excitation. Namely, the edge mode vanishes at 
such k that it has the common tangent with the threshold 
of bulk excitations Eb(k,k y ) with k y = 0. 

The crossings of ID energy bands in Fig. 2 indicates 
the transition from the lowest bulk excitation Eb(k,k y ) 
with k y = to with k y = ± arccos [2mfc/ (A 2 — 4)] . The 
critical momentum k c \ satisfies the condition 



d 2 E h (k cl ,ky)/dk 2 y \ ky= o = 0, (37) 
which gives the solution 

A A 2 

cos fc c i = 1 1 . (38) 

2 4 

The threshold has k y = for k > k c \ > 0. By comparing 
with eq.(35), we find cosfc c i > cosfc m i, which means < 
fc c i < fc rn i. Namely, merging with bulk excitations indeed 
occurs in the range where k y = corresponds to the 
threshold. 

TI-2 regime 

In this case edge modes exist around k = tt. Let us 
assume > when merging occurs at fc = fc m2 > 0. 
The condition for merging is now given by 



-mk 



A 2 



2 + A 

which gives m k = 2 as the solution, or 

cosfc m2 = 1 - A B /2 = 3 - A/2. 



(39) 



(40) 
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zigzag 
edge 




Fig. 3. Zigzag edge lattice ribbon, with two boundaries in (1,1) 
direction. 



In the TI-2 regime, the threshold of bulk excitations oc- 
curs at k y = ir or k y = ± arccos [2mfc/(A 2 — 4)] . It can 
be checked that the bulk energy at k = (/c m 2,7r) be- 
comes the same as the edge mode with the condition (40). 
Hence merging occurs with bulk excitations Eb(k,k y ) 
with k y = 7r. 

The critical momentum k = k C 2 , below which the min- 
imum of Eb(k,k y ) no longer occurs at k y = n, can be 
obtained by condition 



d 2 E h {k c2 ,k y )/dk 2 \ ky 



0, 



(41) 



u 1 Ab 
cos k r ? = 1 — 



(42) 



which gives the solution 

A 2 

2 " 4 

Thus we have the relation cos k c2 < cos k 2m , or k c2 > 
k 2m . Hence merging indeed occurs in the range where 
k y = 7r corresponds to the threshold. In this way, we have 
quantified important characteristics of the edge modes 
shown in Fig. 2. 

4. Zigzag edge 

4-1 Effective one- dimensional form 

Let us now consider a zigzag edge geometry, as illus- 
trated in Fig. 3. Electrons in a zigzag edge geometry are 
confined in a diagonal strip: 1 < y — x < N r , provided 
the edges are placed at y — x = 1 and y — x = N r , normal 
to the (1, — Indirection. The translational invariance re- 
mains along the (1, Indirection, where the conserved mo- 



mentum is given by p = (k x +k y )/y/2. For notational con- 
venience we introduce k = p/y/2 with —n/2 < k < it/2 
and define the new basis set as 
1 



\ K ,J) 



^exp [iK(2I + j)]\I,I + j), 



(43) 



where the site summation goes along the (1,1) direction. 
The phase factor is so chosen that it becomes unity for 
the state | — J, J) . Then, the amplitude $ j in this basis 



satisfies the Schrodingcr equation analogous to eq.(15): 
foj + t u («)*i-i + = Br(«)*j (44) 

where £ has been defined by eq.(7) and hopping matrix 
£ii(k) is given by 



-Ae lK a v + 2 cos fccr 2 



(45) 



We impose the boundary condition: $ = &N r +i = 
0, which is consistent with the zigzag edge geometry. 
Assuming eigenstate of eq.(44) with property $j = 
A$j_i = A 3 $, where |A| < 1, 19 > 20 ) we obtain 



£ + Xt[i(K) + A- i t 11 ( K )J $ = P n (A, k)$ = £ t (k)$ 

(46) 

For later reference purpose, we write the bulk energy 
Eb in terms of variables k = (k x + k y )/2 and £ = (k x — 
k y )/2. From eq.(10) we obtain 

E h (K, £) = ± [2 A 2 (sin 2 k cos 2 £ + cos 2 re sin 2 £) 

1/2 



(As + 4 cos ac cos £)' 



(47) 



4-2 Derivation of spectrum in thermodynamic limit 

We will separate Pn into the Hermitian part Hu and 
the corresponding annihilator Fn. The separation now is 
not straightforward in contrast with the case of straight 
edge. As a preliminary, we introduce the following ma- 
trices: 



ax = (cr x + Vy) /V2, 
cry = (o- y - a x ) /y/2. 



(48) 
(49) 



Then we obtain 
A 

tu (k) = —j= (sin k ax — i cos nay) + 2 cos K a z . (50) 
v2 

We note that the spectrum of each edge mode is an odd 
function of k. Then we introduce a variable 9, which is 
an odd function of k, and make the following transfor- 
mation: 

a$ x = cos 6 ax + sin 8a z , (51) 
og z = cos a z — sin 9 ax , (52) 
cry. They keep the commutation 



and rewrite as ag y 
property: 



in 



2ia ez , 



(53) 



and analogous cyclic ones that are the same as the orig- 
inal Pauli matrices. Then we obtain 

£ = A B (sin 9 a 8x + cos 9 ag z ) , (54) 
\F n + X~Hn = jg x ag x + ie y ag y + ■y Sz ag z , (55) 



where 



.4 



lg x = (A + A ) I — -= sin k cos ( 
V v2 

Jg y = i(X — A^ 1 )— j= COS K, 

v2 



(56) 
(57) 
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lo- 



2 cos k cos 8 



A_ 

71 



sin k sin 



(58) 



We choose $ as eigenstatc of og x since the coefficient 
Ab sin 9 in eq.(54) is an odd function of k. Then terms 
with (JQy and o~g z must combine to form the annihilator. 
Furthermore, since the coefficient of o~g x must be real, 
and be an odd function of k, we require jg x = 0. This 
condition determines 9 in terms of n as 

sm6> = , (59) 



cos# 



^/tan 2 K + 8/A 2 
V&/A 



(60) 

,/tan 2 k + 8/A 2 

H\\ + Fn in eq.(46) as 



In this way, we decompose P\\ 
follows: 

Hn = A B sm8(Tg x , (61) 

F n = {A B COS0 + Jg z )<T 8z +Jg y <Tg y . (62) 

By diagonalizing Hn, the eigenenergy E^(k) is derived 
as 



E t {K) = sA B sin<9, (s = ±l), 



(63) 



where only one of the signs ± is relevant, as derived 
shortly. Note that the spectrum has a form analogous to 
the case of straight edge given by eq.(29). The condition 
for to form the annihilator is given by 



A B cos 9 + j 6z = isjgy 
which determines A for the edge mode as 

1 -A B cos9±y/R 



A±(s) 



with 



R 



2cosk 2/cos6> + s/l/\/2' 



2 cos k (8 cos - 9 — A 



(64) 



(65) 



(66) 



In the case of R < 0, we obtain complex A±(s) with 
absolute value 




0.0 




0.0 




0.0 

Fig. 4. Spectrum of zigzag edge modes. The top panel is for TI-1 
state with A = 1.2, the center is for TI-2 state with A = 4.5. 
The bottom panel shows the boundary case A = 4. All cases 
have A = 1. 



|A±(*)| : 



V$> - sA cos 9 



(67) 



\/8 + sA cos 9 ' 

which is less than unity only for s = 1, and posi- 
tive for A 2 < 8. Therefore, the edge mode must have 
s = 1 in eq.(63), and the group velocity is positive (neg- 
ative) in TP-1 (TP-2) regime. A special case occurs with 
9 = ±7r/2 that corresponds to K = ^1^/2 according to 
eqs.(59) and (60). Actually eq.(65) gives A± = ±i at 
k = tt/2. Thus, at the boundary of the ID BZ, there are 
no edge modes since |A±| = 1. The neighborhood of this 
special point has |A±| < 1, and there should be an edge 
mode. We emphasize that this property is independent 
of A and A, and is specific to the zigzag edge. 

Due to the time-reversal symmetry, we obtain the 
Kramers partner from with the spectrum E^n) = 
— A B sin 9. Figure 4 shows the edge modes together with 
bulk excitations for the TI-1 and TI-2, and the bound- 
ary case A B = 0. The bulk spectrum illustrated is ob- 
tained from Eb(K,l;) as a function of k with fixed £. In 
both regimes, the edge mode becomes gapless at k = 0. 



Hence only k = is the relevant point where a pair of 
edge modes are degenerate by time-reversal invariance. 
We note that the spectrum becomes completely flat with 
A B = 0. 18 > 

In the low momentum region, the spectrum tends to 
the linear dispersion 

Ef^K) = tA^AA b k = t\aA bP . (68) 

Especially, with Ab = ±4, the spectrum becomes the 
same as the corresponding modes in the straight edge. 
Thus we find that the system acquires the axial sym- 
metry only in the case of Ab = ±4 even in the long- 
wavelength limit. This is not surprising since the differ- 
ence in the edge shape remains even for long wavelength. 



4-3 Allowed momentum range for edge modes 

Let us first consider the edge mode in TI-1 regime 
corresponding to right-going edge mode (see Figure 2 
top panel), and A± = A±(l). We restrict to the region of 
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positive k. Since complex A always has a pair of solutions 
with |A±| < 1, the merging momentum K m < n/2 can be 
obtained from the condition 

A+ = 1, (69) 

which is equivalent to 

(8 - A 2 ) cos 2 K m + 2 A B cos n m + A 2 = 0, (70) 

according to eq.(65). Note that there is no solution for 
< cosK m < 1 in the case oi A 2 > 8. The relevant 
solution in the case of A 2 < 8 is given by 

-A B ±^/A 2 B ~A 2 (8 -A 2 ) 



COS t\"rn, 



A 2 



(71) 



The critical value of A B beyond which no real K m exists 
is given by 



-Ay/8 - A 2 < 0, 



(72) 



where we consider only the case < A < 2 as in the the 
straight edge. According to eq.(72), we have three cases: 

(i) no solution for K m < n/2 with As > Ai c ; 

(ii) single n m with A B = A 4c ; 

(hi) two solutions K m \ < K m2 with Ab < A4 C . 

Edge modes in the case of A > 2 will be discussed in §5. 

The threshold of bulk excitations can occur either at 
£ = or £ =/= depending on k. The critical value k c 
separating the two cases is determined by the condition 

d 2 E h (K c ,0/dek=o 
which can be reduced to 



0, 



(73) 



2(A 2 - 4) cos 2 k c - 2A B cos k c — A 
Then we obtain 

A B ± X /A 2 B + 2A 2 (A 2 ~4) 



0. 



COS K r 



(74) 



(75) 



2{A 2 - 4) 

In the case of £ ^ 0, the momentum £ at the threshold 
satisfy the condition 

-2A B cos k 



cos£ = 



(76) 



4+ (4 - A 2 )cos2k' 

At the zone boundary k = n/2, eq.(76) gives £ = ±7r/2 
as the solution. Here the edge mode has the energy 
Ef(n/2) = —A B , and the lowest bulk excitation has the 
same energy E^tt/2, ±7r/2) = —A B - Namely, K = tt/2 
is always a merging point in the zigzag edge for any pa- 
rameter settings. 

In the following, we analyze the spectrum of the edge 
mode near the merging momentum according to classifi- 
cation (i), (ii), (hi) given above. 

Edge modes for the whole BZ 

Let us first consider the case (i): A B > A4 C . With 
A = 1, we obtain A4 C = —y/7, i.e., A c ~ 1.354 from 
eq.(72). Figure 5 shows |A±| and the energy of the edge 
mode relative to the threshold of bulk excitations in this 
case with A = 1.4. 

At the zone boundary, the edge mode merges with bulk 
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Fig. 5. The parameter |A±| and energy difference — _B b be- 
tween edge mode and the threshold of bulk excitations as a 
function of re for A = 1.4 with A = 1. For complex A, wc obtain 
l^+l = l-^-l- I n t ne region of re with real X±, Ef — E\, becomes 
almost zero, but is marginally negative. 



excitations. The difference of energies is expanded as 
E t (n) - E h (K, \) = ^(K - \f + O ((« - \f 



(77) 

Note that only even order terms appear in the expansion 
since both the lowest bulk excitation E\, and edge mode 
energy E^ are symmetric around n/2. 

Edge modes with critical momentum n m 

Next we consider the critical case (ii) characterized by 
single K m with A B ~ A± c . Figure 6 shows |A±| and the 
energy difference. By comparing eqs.(71) and (75), we 
obtain 



< K c l < K m < k c2 < n/2. 



(78) 



Hence we find £ = for the bulk momentum at the 
merging point. 

Edge modes with reentrance 

We finally consider the case (hi): two solutions n m i < 
K m 2 with A B < Ai c . From eqs.(71) and (75) we obtain 
the relation 



< K c l < K m l < K rn2 < k c2 < n/2 



(79) 



Hence we obtain the bulk momentum £ = for both 
merging points. Then we expand E^ — Ei, around merging 
points K m i (i = 1,2) as follows: 

E^(k) - E h (K,0) = «j(k - k„ u ) 2 + O ((k - K m! ;) 3 ) , 

(80) 
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Fig. 6. The same quantities as in Fig. 5 but with A = A c = 
4 — v7 ~ 1.354. Note the scale of the ordinate in the lower 
panel, showing the minute energy difference as compared with 
the overall energy scale of the system. 



Fig. 7. The same quantities as in Fig. 5 but with A = 1. There 
is no edge mode for n between K m \ and re m 2 where one of |A±| 
exceeds unity. 



where the expansion coefficient is calculated as 
A{A i - 8A 2 + A%) 



A 2 A, 



A 2 - 8 ' 



8A 2 + A| 



-, 1/2 



MA 2 



(81) 



Hence we have proven that the threshold of bulk excita- 
tion shares the same energy and velocity with edge mode, 
since the lowest order term of expansion is of second or- 
der. 

At critical value of = A4,,, we obtain K m = K m i = 
K m2 , and both second and third order terms in eq.(80) 
tend to zero. Then the expansion around K m begins from 
fourth order. This explains the nearly flat shape of — 
Eh in Fig. 6. 

Presence of two merging points causes novel phe- 
nomenon in zigzag edge mode. Figure 7 shows the pa- 
rameter |A±| and the energy difference as in previous 
cases. In the region < k < K m \ and K m 2 < K < f, 
we obtain |A±| < 1 and E^ — E h < 0. Namely, the edge 
mode has two separate momentum ranges of its exis- 
tence. The normal edge state (I) starts from zone center 
(k = 0), and vanishes at an intermediate point k = K m i. 
In addition the reentrant part (II) appears near the zone 
boundary: K m 2 < K < \. The reentrant edge mode has 
only marginally lower energy than lowest bulk excitation, 
which approximatively obeys E^ — E^ oc (k — K m 2) ■ In 
the momentum range K m x < K < k,„2, the edge mode 
disappears since one of | A± | exceeds unity. 

Let us summarize the results for different A with fixed 
A shown in Figs. 5, 6 and 7. As A increases from zero, 



the two separate regions for the edge modes widen in 
momentum space simultaneously. Namely, both n m \ and 
7r/2 — K m 2 increase, while n m 2 — n m i decreases with 
increasing A. At critical A c , the two regions merge 
{Kmi = K m 2), and the unbroken edge mode emerges that 
disappears only at the zone boundary. 

4-4 Edge modes in TI-2 regime 

Let us now consider the TI-2 regime with > 0. 
Since the Hamiltonian H(k) has only a linear term of 
Ab, the solution in the TI-2 regime can be obtained from 
that in the TI-1 regime by changing the sign oi Ab- This 
conversion was indeed made in §3 for the straight edge. 

In the zigzag edge, the spectrum given by cq.(63) has 
a negative slope because of As > 0. Corresponding to 
eq.(39), the merging momentum K m > is obtained from 
the condition 



A- = -1, 
instead of eq.(69). Then we obtain 



(82) 



A B ± y/A'j, -A 2 (8- A 2 ) 
cos K m = — ^ — — , (83) 



A 2 



instead of eq.(71), and 



A 4c = Ay? 8 - A 2 > 0, 



(84) 



instead of eq.(72). It is clear that the resultant solution of 
K m is the same as that in TI-1 range with the same |As|. 
Similarly, one can check that all relevant quantities such 
as K m i, "m2, «ci) K c2 also have the same correspondence. 
The expression given by eqs.(77) and (80) remains valid, 
which means that the difference is now positive. This is 
naturally understood since the edge mode E^(k) has the 
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negative slope. 

5. Summary and discussion 

We have analytically obtained spectrum and wave 
function of helical edge states for the BHZ model by 
identifying the annihilator for each case of (1,0) and (1,1) 
edges. The simplicity of the BHZ model has allowed us to 
obtain the complete information of the edge modes. Let 
us finally consider the case A > 2. In the straight edge, 
A± are always real in this case. Except for this differ- 
ence, the property of the edge mode spectrum remains 
the same. In the zigzag edge, the case A > y/8 allows 
no solution for merging momentum n m . This means the 
edge mode is present up to the zone boundary. 

Edge spectrum shows different properties depending 
on edge geometry. In (1,0) -edge case, edge spectrum is 
proportional to sinfc. As the sign of changes from 
negative to positive, the main location of edge mode 
moves from the center of the ID BZ to the boundary. 
This movement is associated with the change of location 
of the bulk energy gap. 

For the (1,1) edge, we have obtained the spectrum in 
the form of ±Ab sin#(ft), where 0(k) is an odd function 
of momentum ft. Since 6(k) does not depend on As as 
shown in eqs.(59) and (60), Ab appears only as the scale 
factor in the spectrum. With Ab = 0, the edge modes 
become completely flat for the whole Brillouin zone, and 
two-fold degenerate as a consequence of the time-reversal 
symmetry. 

The edge mode in the (1,1) geometry contains a novel 
reentrant part with extremely small binding energy. As 
seen from Fig. 7, the binding energy for the reentrant part 
is only 10~ 3 of the overall energy. Except for many-body 
phenomena such as superconductivity and Kondo effect, 
we have been unaware of emergence of such extraordi- 
nary different energy scale. In spite of the tiny binding 
energy, the decay of the wave function toward inside the 
system looks quite normal, as judged by |A±| which de- 
viates clearly from unity. 

Mathematically speaking, the tiny binding energy 
stems from the following factors: 

(i) The zone boundary k = n/2 is a special point where 
the edge mode merges with bulk excitations. 

(ii) Group velocity of the edge mode is the same as that 



of threshold excitation in the bulk at the merging point. 
Let us assume k > for simplicity. The condition (ii) re- 
quires the energy difference — E^ to be proportional to 
(ft — 7r/2) 2 near the zone boundary, and simultaneously to 

[ft ftm2 

) near the critical momentum. Hence in the in- 
tervening region re TO 2 < ft < 7r/2, the growth of E^ — E^ is 
constrained from both ends. Provided 7r/2 — ft m 2 ~ 0.17T, 
we obtain the scaling factor ~ 10~ 2 from the quadratic 
dependence of Ef — E^. It is hoped that more physical 
explanation can be provided in the near future why the 
binding energy is so small. 
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